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Abstract: 

In this article explicit formulas for the recurrence equation 

Pn+l{x) = {AnX + Bn)Pn{x) - CnPn-l{x) 

and the derivative rules 

(t{x)p'^{x) = anPn+l{x) + 0nPn(x) + InPn-lix) 

and 

(j{x)p'n{x) = {an X + Pn) Pn{x) +^nPn-l{x) 

respectively which are valid for the orthogonal polynomial solutions Pn{x) of the differential 
equation 

(j{x) y"{x) + t{x) y'{x) + A„ y{x) = 

of hypergeometric type are developed that depend only on the coefficients (j{x) and t{x) which 
themselves are polynomials w.r.t. x of degrees not larger than 2 and 1, respectively. 

Partial solutions of this problem had been previously published by Tricomi, and recently by 
Yanez, Dehesa and Nikiforov. 

Our formulas yield an algorithm with which it can be decided whether a given holonomic recur- 
rence equation (i.e. one with polynomial coefficients) generates a family of classical orthogonal 
polynomials, and returns the corresponding data (density function, interval) including the stan- 
dardization data in the affirmative case. 

In a similar way, explicit formulas for the coefficients of the recurrence equation and the difference 
rule 

a{x)Vpn{x) = anPn+l{x) + PnPn{x) +JnPn-l{x) 

of the classical orthogonal polynomials of a discrete variable are given that depend only on the 
coefficients a{x) and t{x) of their difference equation 

(j{x) AVy{x) + t{x) Ay{x) + A„ y{x) = . 

Here 

A.y{x) = y(x + 1) — y(x) and Vy(a;) = y{x) — y(x — 1) 

denote the forward and backward difference operators, respectively. In particular this solves the 
corresponding inverse problem to find the classical discrete orthogonal polynomial solutions of a 
given holonomic recurrence equation. 
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1 Polynomials of the Hypergeometric Type 



A long-standing problem in the theory of special functions whose solution can be very helpful 
in applied mathematics as well as in many quantum- mechanical problems of physics [18], 
is the determination of the differentiation formulas of the hypergeometric-type orthogonal 
polynomials Pn{x) only from the coefficients of the differential equation 

a{x) y"{x) + t{x) y'{x) + A„ y{x) = (1) 

which is satisfied by these polynomials 

y{x)=pn{x)=knx'' + ... (nGNo := {0,1,2,... },A;„/0) . (2) 

The coefficients a-{x), t{x) and A„ turn out to be themselves polynomials w.r.t. x of degrees 
not larger than 2, 1 and 0, respectively. 

This problem was partially solved by Tricomi ([21], Chapter IV) in the sense that he was able 
to calculate the coefficients q:„, f3n and 7„ of the derivative rule 

a{x)p^{x) = {an X + ^n) Pn{x) +'ynPn-l{x) ■ (3) 

However, his formula for /3„ was not only in terms of the coefficients of (1) and A;,,,, but 
furthermore the second highest coefficients of Pn{x) were involved, and to evaluate 7„, he 
needed to know also the coefficients of the recurrence equation 

Pn+l{x) = {An X + Bn) Pn{x) - CnPn-l{x) (4) 

another structural property of orthogonal polynomial systems. 

Since the polynomials p„ (a;) given by (2) are completely determined by the differential equation 
and their leading coefficients kn (n G No), it is desirable to obtain the recurrence equation (4) 
and the derivative rule (3) from these informations alone. 

Recently, Yanez, Dehesa and Nikiforov [23] presented such formulas which, however, are ad- 
ditionally in terms of the constant D„, given by a representation of the type 

^''^''> ~ p{x) J {s-x)-+^'^' 
c 

for Pn{x), p{x) being solution of the equation (crp)' = rp, and C being a contour satisfying 
certain boundary conditions. Their development is more general in the sense that they did not 
assume that n is an integer. On the other hand, the assumption that n is an integer implies 
that the contour C is closed, the integral representation (5) being equivalent to the Rodrigues 
represen tat ion 

= Mi) 5^5^ 

where 

n' 

D„ = 5-i;„, (7) 

and the solutions are classical orthogonal polynomials with density p{x). 

In this article, we represent the coefficients of both (3) and (4) in terms of a{x). t{x) and the 

term ratio kn+i/kn alone, hence giving a complete solution of the proposed problem. 
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It is clear that our formulas should depend additionally on the leading coefficients since 
such a standardization can be prescribed arbitrarily. If one takes the monic standardization, 
i.e. kn = 1, then the formulas in fact depend only on the coefficients of the differential equation. 
For the classical orthogonal polynomials our formulas are stronger than Yanez', Dehesa's and 
Nikiforov's result since kn is intrinsic part oi pn{x), whereas the constants Dn,En are not. 
Moreover, we will give Z)„ and En in terms of the coefficients of the differential equation, too. 
Algebraically two identities (differential equation and recurrence equation, e.g.) are needed to 
deduce the third one (derivative rule, e.g.), see [8], whereas here (kind of magic) we would like 
to deduce two from one. That this is possible is due to the analytic knowledge that orthogonal 
polynomial solutions of the differential equation (1) satisfy some structural properties, namely, 
the recurrence equation and derivative rule take special forms. 

We make the general assumption that our polynomials Pn{x) are orthogonal w.r.t. a measure 
iu, i.e. 

if m / n 

Pn{x) Pm{x) dnix) = { (8) 



/ 



hn 7^ if m = n 



where / denotes an appropriate integration path, for example a real interval. 

Major tools in our development are the following well-known structural properties of such 

families of orthogonal polynomials. 

Lemma 1 Any system of polynomials {pnix) \n ENq}, Pn being of exact degree n, orthogo- 
nal with respect to a measure n, satisfies a three-term recurrence equation of the form (4) 

Pn+l{x) = {An X + Bn) Pn{x) - CnPn-l{x) (n G Nq , P-l{x) = 0) , 

An,Bn and C„ not depending on x. 

Proof: This property is well-known (see e.g. [21], Chapter IV). To prove it, one substitutes 
(2), equates the coefficients of x""*"^, and gets immediately that 

An = ^. (9) 

With this choice, we study the difference Pn+i{x) — AnXpn{x). Since this is a polynomial of 
degree not larger than n, it can be decomposed as 

n 

Pn+l{x) -AnXPn{x) = ^djPn{x) . 

j=0 

We choose m < n — 2 and multiply by pm{x). Integrating with respect to /j, yields 

where on the right hand side (8) was applied. Both left hand integrals vanish since Pm{x) is 
orthogonal to Pn+i{x) and since xpm{x) as a polynomial of degree not larger than n — 1 is 
orthogonal to Pn{x), implying dm = 0. This gives the result. □ 

The second important structural property for our considerations is given by 
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Lemma 2 Any system of polynomials {pn{x) | n E No}, Pn being of exact degree n, that 
are solutions of the differential equations (1) and furthermore orthogonal with respect to a 
measure ij,{x) = p{x) dx having weight function p[x) > 0, satisfies a derivative rule of the 
form (3) 

o{x)p'j^{x) = {anX + (5n)Pn{x) + %Pn-i{x) (n G N := {1,2,3,...}) , 
an,Pn and 7„ not depending on x. 

Proof: Substituting (2) and equating the coefficients of x""*"^, one gets immediately that 

ctn = an . (10) 

In [13], § 5 it is shown by an elementary argument that under the given conditions the solutions 
Pn{x) of the differential equations (1) are orthogonal with respect to the weight function 

p(a;) := -^e^^'^'' > , (11) 
a{x} 

given by Pearson 's differential equation 

■^(^a{x) p{x)^ = t{x)p{x) 

for a suitable constant C, in a suitable interval I (depending on the zeros of cr{x)). Hence 
multiplying (1) by p(x), the differential equation takes the selfadjoint form 

^ (a{x) p{x) y'{x)^ + A„ p{x) y{x) = . 

Using this identity, Tricomi showed that ([21], IV (4.10)) 

I a{x)p{x)p'^{x)f{x)dx = (12) 

for any polynomial f{x) of degree < n — 2. If (10) holds, then the degree of a{x) p'^ix) — cunX 
is < n. Hence one can write 

n 

a{x)p'n{x) -a.nX = ^ejPn{x) . 

As above, from (12) one can deduce that ej = for < j < n — 2 (see [21], Chapter IV). □ 
An immediate consequence is the following 

Corollary 1 Any system of polynomials {pn{x) | n G Nq}, being of exact degree n, that 
are solutions of the differential equation (1) and furthermore orthogonal with respect to a 
measure p{x) = p{x) dx having weight function p{x) > 0, satisfies a derivative rule of the form 

a{x)p'^{x) = anPn+i{x) + (3nPn{x) +7„p„-i(a;) (n G N) , (13) 

an,Pn and 7„ not depending on x. 
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Proof: Substituting (2) in (13), and equating the coefficients of x""*"^, one gets immediately 
that ^ 

an = an—^ . (14) 

Substituting (4) in (3) one gets moreover 

a{x)p^{x) = {anX + Pn)Pnix) +'ynPn-lix) 

^ (^Pn+l{x) - BnPn{x) + CnPn-l{x)^ + Pn{x) +'^nPn-l{x) , 

hence (13) is vahd with 

— a — a _ ~ — ~ I ~ i-i 

Oln — A : Pn — Pn OLn , i In — Tn ~r , ■ I— ' 

■^n -^n -^n 

2 Classical Orthogonal Polynomials of an Interval 

In this section we give the proposed expUcit recurrence equation and derivative rule formulas. 
Assume a family of differential equations (1) is given for n G No, with continuous functions 

a{x),T{x), and constants A„, and we search for polynomial solutions (2) of degree n. Then 
since pi(x) is linear, one deduces that t{x) must be an at most linear polynomial, and since 
P2{x) is quadratic, one deduces that a{x) must be an at most quadratic polynomial [3]. Hence 
we may assume that 

(j{x) := ax^ + bx + c , t{x) := dx + e . (15) 

Equating coefficients of the highest powers in (1) for generic Pn{x), given by (2), one 
deduces that moreover 

an{n — 1) + c?n + A„ = or A„ = —{an{n — 1) + dn) . (16) 

Hence only if the differential equation takes the special form 

{ax^ + bx + c) y"{x) + {dx + e) y'{x) - {an{n - 1) + dn) y{x) = , (17) 

it can have polynomial solutions. 

Moreover we can assume that A„, / for n G N, hence a{n — 1) + d 7^ for n G N since 
otherwise no orthogonal polynomial solutions can exist. This is discussed in detail in [13]. In 
particular, d 7^ 0. 

In the following theorem, we give explicit representations of the corresponding recurrence 
equation and derivative rule in terms of the given a, 6, c, c?, e and the term ratio A;„_|_i/A;„. 

Theorem 1 Let Pn{x) = A;„ + . . . (n G No) be a family of polynomial solutions of the 
system of differential equations (17) that are orthogonal with respect to a weight function 
p{x). Then the derivative rule (13) 

a{x)p'n{x) = anPn+i{x) + (3nPn{x) +7„p„-i(a;) 

is valid with 

an= an 
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_ n{a{n - 1) + d) {bd - 2ae) 
^" ~ (2a(n-l) + d) (2an + d) ^ ^ 

n (a(n— 1) + d)(a(n — 2) + (i)(n (an + d) (4ac — 6^)+ae^+cc?^— 6c?e) A,'„ 
" (a(2n-l)+d) (a(2n-3)+d) (2a(n-l) + d)2 ' ^ ^ 

and the recurrence equation (4) 

Pn+lix) = {AnX + Bn)Pnix) - CnPn-lix) 



is valid with 



n 



and 



7„ being given by (19). 



^ _ kn+i 2bn {a{n - 1) + d) + e{d - 2a) 

kn ' (2a(n - I) + d) {2an + d) ^ ' 

" kn a{n — 1) + d ' 



Proof: The values of An and were already obtained in Lemma 1 and Corollary 1 . 
By Lemma 1 the polynomials satisfy a recurrence equation of type (4): 

Pn+l{x) = {An X + Bn) Pn{x) - CnPn-l{x) ■ (22) 

Next, we differentiate (22) twice and get 

P'n+lix) = AnPnix) + {A^ X + B^) p'^ix) - Cnp'n-lix) (23) 

and 

p';,+^ix) = 2Anp'nix) + iAnX + Bn)p';,ix) - Cnp'n-lix) ■ 

We multiply the last equation by a{x) 

(j{x)Pn+i{x) = 2An(j{x)p'n{x) + {AnX + Bn)'j{x)Pn{x) - Cn (rix) p'^,_i{x) 

and use the differential equation to replace the second derivatives by those of lower order 

-(t{x)p'^_^_^{x) + Xn+l Pn+l{x)^ = 2An (t{x) p'^ix) - {An X + Bn) (t{x)p'^{x) +A„p„(x)j 

+ Cn (T{x)p'n_i{x) + A 

After substituting (23) on the left hand side, and subtracting r {Ax + B)p'n — rCnp'n-i^ 
arrive at 

-t{x) AnPn{x) - \n+lPn+l{x) = 2An a{x)p'^{x) - {AnX + Bn) KPn{x) + C„ A„_ip„_i(2;) . 

Next, on the right hand side, we replace {An x + Bn)pn{x) by pn-\-i{x) + CnPn-i{x) according 
to (22), and get 

-t{x) AnPn{x) - Xn+lPn+l{x) = 2An a{x) p'n{x) - XnPn+l{x) + (A„_i - A„) CnPn-l{x) , 
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or rewritten 

l)Pn+l{x) = An {2a{x)p'n{x) +T{x)pn{x)) + (A„_i 

~ A„) Cn Pn-l{x) . 

Now we substitute the representation of Corollary 1 

a{x)p'n{x) = anPn+l{x) + (3nPn{x) +'ynPn-l{x) 

to deduce 

^(A„ - A„+i) -2 An On^ Pn+l{x) = An (2 (3n+T{x)) Pn{x) + {2 A^ 7„ + (A„_i - A„) C„) Pn-l{x) 

after subtracting 2y4„ a„p„+i(a;). Replacing pn-\-i{x) according to (22), we arrive at the 
identity 

[An {2l3n + T{x)) + (2^„a„ - (A„ - A„+i))(^„ a; + p„(a;) = 

-(((An - A„+i) -2Anan) C„ + 2^„7„ + (A„_i - A„) Cn)pn-i{x) . 

Since Pn{x) is a polynomial of exact degree n, this relation can only be valid if the coefficients 
of both Pn{x) and Pn-i{x) vanish, since otherwise the polynomial on the left hand side has 
degree > n whereas the polynomial on the right hand side has degree n — 1 , a contradiction. 
The coefficients must vanish as polynomials in a;, and equating coefficients we are led to the 
three equations 

An (A„+i -\n+d + 2an An) = . 

Ane + 2AnBnan + 2Anfin — Bn A„ + En A„+i = 

and 

2 An CnOtn — 2 An 7n — C„ A„_i + Cn A„_|-i = . 

Whereas the first of these equation does not contain any news but restates a relationship 
between y4„, Q!„ and A„, the second and third of these equations (using (9), (14) and (16)) 
can be rewritten as 

= d K' ^ ^ 

and as (21). Hence S„ and C„ are known as soon as /3„ and 7„ are. 

We finally need two more equations to find /3„ and 7„. To deduce one of these equations, and 
to find fin, and hence S„, we substitute 

Pn(*^) — X ~\~ kn X -\- kn X -\- . . . (^^) 

in the three equations considered, namely the differential equation, the recurrence equation 
and the derivative rule. As we already saw, equating the coefficients of the highest powers of 
X yields (16), (9) and (14). If we equate the coefficients of the next highest powers of x, we 
get three more equations, involving two more variables though, namely k'n and These 
are the equations 



bnkn — enkn — bn^ kn — 2akn + dk'n + 2ank'n = , (26) 



Bn kn + An kn ^n+l ~ ^ 5 (2'^) 



and 

bnkn — (3nkn — ak'n + ank'n — CXn k'n^i = . (28) 
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Equation (26) immediately gives 



^ n {b{7i - 1) + e) 
kn 2a{n-l)+d ' ^ ' 

whereas from (27)-(28) one can eliminate /c^+i- This gives a second equation between Bn and 
/3n which together with (24) and (29) yields (18) and (20). 

To deduce 7^1, we equate the coefficients of the next highest powers in the differential equation, 
recurrence equation and derivative rule, introducing two more auxiliary variables k!^ and A;^_|_i 
which can be eliminated. This procedure generates one more equation between C„ and 7„ 
finally deducing (19). □ 

Note that the results given in Theorem 1 can also be deduced completely automatically by 
elimination methods based on Grobner basis calculations. With the computer algebra systems 
Maple and REDUCE we were successful doing so. For the purpose of finding C„, we 

substitute (25) in the differential equations for Pn{x) and for p„_|_i(a;), and in the recurrence 
equation. Equating the three highest coefficients in any of these three equations yields nine 
nonlinear equations in the nine unknowns 

-^m Bji, Cjn Ajj, Xn+lj ■ 

By a Grobner basis computation (invoked by the solve command of the utilized computer 
algebra system) it turns out that there is a unique solution, given by Theorem 1, see also 
Corollary 4 and (33). Note that therefore the formulas for An,Bn, and C„ of Theorem 1 are 
valid without the hypothesis of a weight function p{x). 

Similarly, to find Q-n-, l^niln^ we substitute (25) in the differential equations for Pn{x) and for 
Pn+i{x)^ and in the derivative rule. Equating the three highest coefficients in any of these 
three equations yields nine nonlinear equations in the nine unknowns 

(^m Pn: Tni An, A^+i, A;^^, A;„_|_]^, kji, 5 

and a Grobner basis computation generates the unique solution, given by Theorem 1. Note 
that we were not able to separate the two problems in a similar way based on hand calculations. 

Our theorem has immediate consequences. 

Corollary 2 Let Pn{x) = knX^ + ... (n G No) be a family of polynomial solutions of the 
system of differential equations (17) that are orthogonal with respect to a weight function 
p{x). Then the derivative rule (3) 

a{x)p'^{x) = {an X + ^n) Pn{x) +%Pn-l{x) 

is valid with 

(in = an 
~ _ {ab{n-l) -ae + bd)n 

2a(n-l) + d ^^^^ 

^ -a c + ^sn 

a [n — 1) + a 

In and Cn being given by Theorem 1. 
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Proof: Substituting (4) in (13) yields (3) with 

This yields the result. □ 

Note that Theorem 1 describes the variety of different recurrence equation formulas known 
in the literature ([1], 22.7) by one single formula. Similarly all the different derivative rule 
formulas ([1], 22.8) are governed by a single formula through Corollary 2. 
Theorem 1 shows in particular that the recurrence equation and derivative rule can be obtained 
by purely rational arithmetic whenever 

G Q(n) , 



i.e., if kn is a hypergeometric term. This is obviously true if kn = 1, i.e., in the monic case. 
But also all other standardizations that are used in practice (see e.g. [1], Chapter 22) are of 
this type.^ 

In the case of the orthonormal standardization given by 

hn = 1 

it is not in general true that kn is a hypergeometric term. On the other hand, if kn is a 
hypergeometric term, hn inherits this property. 

Corollary 3 Let Pn{x) = A;„ + . . . (n G No) be a family of orthogonal polynomial solutions 
of the system of differential equations (17). Then the relation 

hn+i {n + 1) {an + d){a{n — l) + d) f 6(n+l) + e/ 7 ^' 



hn (a(2n + 3)+c?)(a(2n + l)+c?) 
is valid. 

Proof Tricomi ([21], IV (2.2), see also [1], (22.1.5)) proved that 

n 



An application of Theorem 1 yields (32). □ 

Next we would like to give a general formula for the term ratio of the coefficients k'^ in terms 
of the given term ratio of kn. 

Corollary 4 Let Pnix) = knx'^ + k\^ x'^^^ + . . . (n G No ) be a family of orthogonal polynomial 
solutions of the system of differential equations (17). Then the relation 

k'n+i _ n + 1 {hn + e)(2a(n - 1) + d) kn+i 
k'n n {h{n — I) + e){2an + d) kn 

is valid. 



^Only in one instance, this is not so: For the Chebyshev polynomials T„(a;) one has kn+i/kn = 2 (n G 
and ki/ko = 1. If one redefines To{x) := 1/2, then kn+i/kn =26 Q(n). 
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Proof: This follows from 



kL kn4-^ kn k' 



using (29). □ 
Similarly, one gets 

k'^ n{n-l){n'^b'^ -3nb'^ + 2nbe + 2cna-2ca-3be + 2b'^ + 6"^ + cd) 



2 (2 an - 2 a + d) {d - 3 a + 2 an) 



(33) 



also deduced by the automatic elimination method mentioned before, and a similar equation 
for k'^j^ijk'^, see the Appendix. 

We furthermore obtain the term ratio of the numbers of an integral representation of 
type (5) considered in [23]. Note that in the orthogonal polynomial case the contour C is 
closed, and hence by Cauchy's integral formula representation (5) is equivalent to a Rodrigues 
representation (6) with 



Dn — TT^ En 
2m 



We get for and D„, respectively 



Corollary 5 Let Pn{x) = kn x'^ + . . . (n G No) be a family of orthogonal polynomial solutions 
of the system of differential equations (17). Then Pn{x) have a Rodrigues representation (6) 
and an integral representation (5) with closed contour C surrounding s = x, and one has for 



Efi and D„ the relations 



Dn+i (n + 1) (a(n - 1) + d) k 



'n+l 



Dn (a(2n - 1) + d) {2an + d) 

and 

En+i {a{n -l)+d) 



En {a{2n - I) + d){2an + d) k, 
Proof: In ([23], (13)) it was shown that 

1 Dn (n + 1) {a{n - 1) + d) 



An Dn+i (a(2n - I) + d) {2an + d) ' 

An application of (9) leads to the term ratio for D„. The term ratio for En follows then from 
(7). □ 

Note that Corollary 5 again describes all the different Rodrigues formulas ([1], 22.11) known 
in the literature by one single formula. 

It is well-known ([3], see also [4], [13]) that polynomial solutions of (1) can be classified accord- 
ing to the zeros of (j{x), leading to the normal forms of Table 1 besides linear transformations 
X 1-^ Ax + B. The type of differential equation that we consider is invariant under such a 
transformation. Orthogonal polynomial solutions according to this classification exist if and 
only if the function 

C f^dx 

a{x} 

given by (11) yields a weight function in the interval given by the zeros of (j{x), i.e. the 
corresponding integrals converge and p{x) > for some C. 
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1. a = b = c = e = 0, d = 1 =^ Pnix) = a;", 

2. a = b = e = 0, c = 1, d = —2 =^ Pn{x) = Hn{x), the Hermite polynomials, 

3. a = c = 0, &=1, c?=— 1, e = Q:+l = L^"^ (x), the Laguerre poiynomiais, 
4a. a = l, 6 = c = c? = e = 0, = x", 

4b. a = l, 6 = c = 0, d = a + 2, e = 2 =^ Pn{x) = Bi^\x), the Bessel polynomials, 

5. a = l,b = 0,c = — l,d = a+f3+2,e = a—f3 =^ Pn{x) = Pn"''^\x), the Jacobi polynomials. 
Table 1: Normal Forms of Polynomial Solutions 



This shows that the only orthogonal polynomial solutions are linear transforms of the Hermite, 
Laguerre, and Jacobi polynomials, hence using a mathematical dictionary one can always 
deduce the recurrence equation and derivative rules. Note, however, that this approach (in 
general) requires the work with radicals, namely the zeros of the quadratic polynomial (j[x), 
whereas our approach is completely rational: Given kn+i/kn G Q(n), the recurrence equation 
and derivative rules are given rationally in Theorem 1. 

Note that the formulas of Theorem 1 are also valid for the Bessel polynomials ([18], p. 24) 



Bi-\x) 



(2n)!a^ 
n!2" 



1^1 



-n 
-2n 



g2/x 



2" dx 



This is so since the Bessel polynomials do satisfy both a recurrence equation and a derivative 
rule of the desired type (see e.g. [23]), despite the fact that the corresponding function 



p{x) 



C f^dx ^ a -2lx 



o{x) 



does not constitute a weight function on the real axis. The validity of both a recurrence 
equation and a derivative rule of the given types, however, was the only assumption in the 
proof of Theorem 1. 

Although the Jacobi polynomials Pi°''^\x) do only constitute orthogonal polynomials for 
a, /3 > — 1, by a simple argument it can be shown that the structural properties like recurrence 
equation and derivative rule remain valid for arbitrary values of a, [3. A similar comment 
applies to the other parameterized families of Table 1. Hence Theorem 1 is valid also in these 
cases. 

Theorem 1 is even valid in the case of Table 1:4a, and its recurrence equation part also for 



Table 1:1 with the trivial solution 



In both cases we have the recurrence equation 



Pn+i{x) = xpn{x), and in the first case we receive the derivative rule x p'^{x) = np„_|_i(a;). 
Note that there is another derivative rule xp'^{x) = npn{x) which cannot be discovered by 
Theorem 1. 



In the next section we will use the fact that these equations are given explicitly to solve an 
inverse problem. 
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3 The Inverse Characterization Problem 



Assume you have a polynomial system given by a differential equation (1). Then by the 
classification of Table 1 it is easy to identify the system. On the other hand, given an arbitrary 
holonomic three-term recurrence equation 

qn{x) Pn+2{x) + rn{x) Pn+i{x) + Sn{x) Pn{x) = (g„(a;), r„(a;), G Q[n, x]) , (34) 

it is less obvious to find out whether there is a polynomial system 

P„(a;) = A;„a;" + ... (nGNo,A;„/0) 

satisfying (34), being a linear transform of one of the classical systems (Hermite, Laguerre, 
Jacobi, Bessel), and to identify the system in the affirmative case. In this section we present an 
algorithm for this purpose. Note that Koornwinder and Swarttouw have also considered this 
question and propose a solution based on the careful ad hoc analysis of the input polynomials 
g„,r„, and s„. Their Maple implementation works for a part of the so-called Askey- Wilson 
scheme ([2], see also [12]). 

Let us start with a recurrence equation of type (34). We assume that neither qn-i{x) nor 
Sn{x) has a nonnegative integer zero since otherwise this recurrence equation cannot be used 
to determine Pn{x) iteratively from Po{x) (with P-i{x) = 0) for all n > 1 or is worthless in 
the backward direction. Define 

{0, if neither qn-i{x) nor Sn{x) have a nonnegative integer zero 
max{n G No | n is a zero of either qn-i{x) or Sn{x)} + 1, otherwise 

Then we consider Pn{x) := Pn+Nix) instead of Pn{x), see § 4 for an example with N > 0. In 
this situation we rewrite (34) by substituting n hj n + N and replacing Pn{x) by Pn-N{x)- 
For simplicity we rename g„(a:;), r„(a;) and Sn{x)^ and assume in the sequel that the recurrence 
equation 

qn{x)pn+2{x) + rn{x)pn+i{x) + Sn{x)pn{x) =0 (g„ (s) , r„ (x) , s„ (s) GQ[n,x]) (35) 

is valid, but now neither qn-i{x) nor Sfii^x^ have nonnegative integer zeros. We search for 
solutions 

Pr,{x)=Kx'' + ... (nGNo,A;„/0) (36) 
which reads in terms of the original family P„ (n = TV, TV + 1, . . .) 

Pn+N{x) =knX^ + ... (wGNq). 

Next, we divide (35) by qn{x), and replace n by n — 1. This brings (35) into the form 

Pn+l{x) = tn{x) Pn{x) + Unix) Pn-l{x) {tn{x),Un{x) & Q{n,x)) . (37) 

For Pn{x) being a linear transform of a classical orthogonal system, there is a recurrence 
equation (4) 

Pn+l{x) = {AnX + Bn)pn{x) -CnPn-l{x) (^„, S„, C„ G Q(n), ^„ / 0) , (38) 
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therefore (37) and (38) must agree. We would like to conclude that tn{x) = AnX + Bn-, and 
Unix) = —Cn which follows if we can show that Pn{x)/pn-i{x) ^ Q{n,x). To prove this 
assertion, we assume that 

-^^,eQ{n,x). 

Pn-l{x) 

Hence there are P{n,x) G Q[n][a;] and Q{n,x) G Q[n][a;] with gcd^ {P{n,x),Q{n,x)) = 1 such 
that the relation 

Q{n,x)pn{x) = P(n,x)pn-i{x) (39) 

holds. It is a classical result for orthogonal polynomials that gcdr^{pn{x),Pn-i{x)) = 1 since 
their zeros separate each other (see e.g. [21], IV. 6). Hence from (39) we conclude that 

P(n,x) = SnPn{x) and Q{n,x) = SnPn-i{x) . 

Since by assumption P(n, x) G Q[n][a;] should be a polynomial of fixed degree with respect to 
X, and since Pn{x) has degree n, this gives an obvious contradiction. 

Therefore we can conclude that tn{x) = An x + Bn, and Un{x) = —Cn- Hence if (37) does not 
have this form, i.e., if either tn{x) is not linear in x or Un{x) is not a constant with respect to 
X, we see that Pn{x) cannot be a linear transform of a classical orthogonal polynomial system. 
In the positive case, we can assume the form (38). 

Since we propose solutions (36), equating the coefficients of a;""*"^ in (38) we get 

^=Ar, = ^ {Vn,Wr,eQ[n]) . (40) 

Kn Wn 

Hence the given A^ = v^/wn G Q(?t.) generates the term ratio kn+i/kn, and in particular 
kn turns out to be a hypergeometric term which is uniquely determined by (40) up to a 
normalization constant = po{x). Since the zeros of Wn correspond to the zeros of qn-i{x), 
kn is defined by (40) for all n G N from ko. 

In the next step we can eliminate the dependency of kn by generating a recurrence equation 
for the corresponding monic polynomials Pn{x) = Pn{x)/kn- For Pn{x) we get by (40) 

Pn+l{x) = {x + ^\ Pn{x) - Pn-l{x) = (x + Bn) Pn{x) - CnPn-l{x) 

\ An J AnAn-l ' 

with 

Bn = ^e Q(n) and (?„ = ^" G Q(n) . 

Then our formulas (20)-(21) read in terms of Bn and Cn 

~ ^ 2bn {n{n - 1) + d) + cjd - 2a) 

{2a{n - I) + d) {2an + d) ^ ' 

and 

(a(2n-l)+d)(a(2n-3)+d) V (2a (n-l)+d)2 V > ^ >))' ^ > 

and these are independent of kn by construction. 

Now we would like to deduce a,b,c,d and e from (41)-(42). Note that as soon as we have 
found these five values, we can apply a linear transform (according to the zeros of cr{x)) to 
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bring the differential equation in one of the forms of Table 1 which finally gives us the desired 
information. 

We can assume that S„ and Cn are in lowest terms. If the degree of either the numerator or 
the denominator of ^„ is larger than 2, then by (41) Pn{x) is not a classical system. Similarly, 
if the degree of either the numerator or the denominator of Cn is larger than 4, by (42) the 
same conclusion follows. 

Otherwise we can multiply (41) and (42) by their common denominators, and bring them 
therefore in polynomial form. Both resulting equations must be polynomial identitites in 
the variable n, hence all of their coefficients must vanish. This gives a nonlinear system of 
equations for the unknowns a, b, c, d and e. Any solution of this system with not both a and d 
being zero yields a differential equation (17), and hence given such a solution one can decide 
whether the corresponding solutions Pn{x) are generated by a density (11). Therefore our 
question can be resolved in this case. 

If the nonlinear system does not have such a solution, we deduce that no such values a, 6, c, d 
and e exist, hence no such differential equation is satisfied by Pn{x), implying that the system 
is not a linear transformation of a classical orthogonal polynomial system. 
Hence the whole question boils down to decide whether the given nonlinear system has non- 
trivial solutions, and to find these solutions in the affirmative case. As a matter of fact, with 
Grobner bases methods, this question can be decided algorithmically [15]-[17]. Such an algo- 
rithm is implemented, e.g., in the computer algebra system REDUCE [16], and Maple's solve 
command can also solve such a system. 

Note that the solution of the nonlinear system is not necessarily unique. For example, the 
Chebyshev polynomials of the first and second kind Tn{x) and Un{x) satisfy the same recur- 
rence equation, but a different differential equation. We will consider this example in more 
detail later. 

If we apply this algorithm to the recurrence equa.tion pn+2{x) — x Pn+i{x) of the power p„ (a;) = 
a;", it generates the complete solution set, given by Table 1:1 and 1:4a. 
The following statement summarizes the above considerations. 

Algorithm 1 This algorithm decides whether a given holonomic three-term recurrence equa- 
tion has classical orthogonal polynomial solutions, and returns their data if applicable. The 
algorithm is applicable to all entries of Table 1 independently of the orthogonality of the 
system under consideration. 

1. Input: a holonomic three-term recurrence equation 

Qn{x) Pn+2{x) + rn{x) Pn+l{x) + Sn{x) Pn{x) =0 {q 

n {x^ , r^i (s) , Sjj {x) G Q[n,x]) . 

2. Shift: Shift by max{n G Nq |n is zero of either qn-i{x) or s„(a;)} + 1 if necessary. 

3. Rewriting: Rewrite the recurrence equation in the form 

p„+i{x) = tn{x)pn{x) + Un{x) pn-i{x) (i„ (s) , it„ (s) GQ(n,a;)) . 

If either tn{x) is not a polynomial of degree one in x or Un{x) is not constant with respect 
to then return "no classical orthogonal polynomial solution exists"; exit. 

4. Standardization: Given now A,,, i?„ and C„ by 

Pn+l{x) = {AnX + Bn)Pn{x) - CnPn-l{x) (^„ , S„ , C„ G Q(n) , ^„ / 0) , 
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define 

according to (40). 

5. Make monic: Set 

Bn:=^^ Q(n) and (?„ := ^" G Q(n) 

and bring them in lowest terms. If the degree of either the numerator or the denominator 
of Bji is larger than 2, or if the degree of either the numerator or the denominator of Cn 
is larger than 4, return "no classical orthogonal polynomial solution exists"; 
exit. 

6. Polynomial Identities: Set 

~ _2bn{a{n-l)+d)+e{d-2a) 
" ~ (2a(n -l)+d) (2an + d) 

and 

~ -niain-2)+d) ( 6(n-l) + e /, , ,,\\ 

Cn = -1—, , ; ;r c+ -, — -, -77 (ae-M) -flfe (n- 1) , 

(a(2n-l)+d)(a(2n-3)+d) V (2a (n-l)+d)2 V > ^ >))' 

using the unknowns a, 6, c, d and e. Multiply these identities by their common denomi- 
nators, and bring them therefore in polynomial form. 

7. Equating Coefficients: Equate the coefficients of the powers of n in the two resulting 
equations. This results in a nonlinear system in the unknowns a, 6, c, d and e. Solve 
this system by Grobner bases methods. If the system has no solution, then return 
"no classical orthogonal polynomial solution exists"; exit. 

8. Output: Return the classical orthogonal polynomial solutions of the differential equa- 
tions (17) given by the solution vectors (a, 6, c, rf, e) of the last step, according to the 
classification of Table 1, together with the information about the standardization given 
by (40) . This information includes the density 

p{x) 1 f ^ dx 
C a{x) 

given by (11), and the interval by the zeros of a{x). □ 



We would like to give the following comments on the above algorithm: 

1. The use of Grobner bases is not always necessary. The following observation yields an 
ad hoc method to solve the nonlinear system. Observe that the coefficients of the powers 
of n of the polynomial identity concerning S„ of step 5 of the algorithm can be written 
using the variables 

{a'^,ad,ae,ab,d'^,de,db} (43) 
Then all the derived equations are linear in the seven variables (43). 
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Furthermore the coefficients of the powers of n of the polynomial identity concerning C„ 
of step 5 of the algorithm essentially are products of exactly four terms out of a, b, c, d, e. 
Any of these can be written as a product of two of the variables 

{di = a^,d2 = ad,d3 = ae,d4^ = ab,d5 = d^,de=de,dY = db,ds = ac,dg = dc} . (44) 

This is the set of variables (43) plus the two variables ds,dQ. All these equations are of 
second order in the nine variables (44). 

Obviously the resulting system can be solved by first finding the solution space corre- 
sponding to the linear subsystem, which then can be substituted in the second order 
subsystem. The resulting second order system can be solved by ad hoc elimination (and 
possibly rational factorization) . 

If one has found the variables given by (44), then it is easy to calculate a,b,c,d and e, 
or one realizes that no such solution exists. 

2. Note moreover that, although Grobner bases techniques apply only rational arithmetic, 
hence give rational solutions only, the technique described shows that solutions are also 
detected if they involve radicals. 

3. If the input recurrence equation has further parameters, in step 6 of the algorithm one 
should solve for all variables including these additional ones, see Example 2. 

Example 1 As a first example, we consider the recurrence equation 

{n + 2)Pn+2{x)-x{n + l)Pn+i{x)+nPnix) = . 

Since so{x) = 0, we see that the shift Pn{x) := Pn+i{x) is necessary. For Pn{x), we have the 
recurrence equation 

(n + 3)pn+2{x) -x{n + 2)pn+i{x) + {n + l)pn{x) = . (45) 

In the first steps this recurrence equation is brought into the form 

n + 1 n 

P„+l{x} = —XPn{x} -^Pn-l{x) , 

n + 2 n + 2 



hence 



and therefore 



_ kn+l _ n+l _ Vn 
kn n + 2 Wr, 



kn — ^ ko 



n + 1 

Moreover, for monicp„(a;) =Pn{x)/kn we get 

P„+l{x) = XPn{x) +Pn^i{x) , 

hence Bn = and C„ = 1. In step 5 of the algorithm, the polynomial identity concerning B„ 
then reads as 

-2abn^ + {2ab-2db)n-de + 2ae = , 
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leading to the linear system 

ab = 0, 2ab-2db = 0, de-2ae = 

for the variables ab,db,de,ae, with the solution 

{de = 2ae,ab = 0,db = 0} , (46) 

ae being arbitrary. After substituting the corresponding equations 

{4 = 2c?3,c?4 = 0,c?7 = 0} (47) 

into the polynomial identity concerning Cn, we equate the coefficients and receive the second 
order equations 

4di (4di + d8) = , (48) 

23 di 4 - 28 di d2 + 12 di"^ - 8d2 ds + d^'^ = , (49) 

92 c?i^ - 96 di d2 + 24 4 + 5 ^2 dg - 20 di dg + 20 di dg + rfs^ = , (50) 

92 di c?2 - 56 c?i^ - 48 c?i c?5 + 8 d2 4 - 6 ^2 dg + ds dg + 12di dg - 8 di dg = , (51) 

-8di (8^1-4^2 -f^9 + 2^8) =0 (52) 

in terms of the variables (44). The first of these equations leads to two possibilities: either 
c?i = or c?8 = — 4c?i. One realizes quickly that the first of these possibilities implies a = d = 
which is not allowed. Hence we must have 

c?8 = — 4c?i, or c = — 4a, (53) 

and di ^ 0, i.e. a ^ 0. At this point we have already determined a{x) since by (46) one has 
a 6 = 0, hence 6 = and therefore 

a{x) = ax^ + bx + c = a{x^ — 4) . 

Hence possible orthogonal polynomial solutions of (45) are defined in the interval [—2,2]. 
We substitute now (53) in (49)-(52). Then the last equation reads as 

-8c?i (-4c?2 - dg) = . 

Since c?i / 0, we conclude that 

dg = -4d2 or c = -4 a . (54) 

In terms of a, b, c, d and e this yields nothing new, but it shows the compatibility of (52) with 
(48). 

Substituting (54) in (48)-(52) gives two trivial identities, and three complicated ones. In these 
three equations, we finally resubstitute the original variables by (44), and after a rational 
factorization we get 

{Sa-d){a-d){2a-df = 0, 
—4 a [3 a — d) {a — d) [2 a — d) = 0, 
(|l2a2 _ 16ad + 4d2 + 6^) = . 
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Hence either 

d = a, or d = 2a, or d = 3a . 
In the first of these cases, one gets e = and the differential equation 

(x^ - 4) y"{x) + X y'{x) - n{n - 2) y{x) = (55) 

corresponding to the density 

1 r 44 1 



1 r44dx 

p{x) = -— eJ ''(^) 



The corresponding orthogonal polynomials are multiples of translated Chebyshev polynomials 
of the first kind 

Pn{x) = K Cn{x) = Cn{x) = ^ T„ (a;/2) (n > 0) (56) 

n + 1 n + 1 

(see e.g. [1], Table 22.2, and (22.5.11); Cn{x) are monic, but Co = 2, see also Table 22.7), 
hence finally 

Pn{x) =Pn-l{x) = ^Tn-l{x/2) (n > 1) . 

n 

In the second of the above cases, i.e. for d = 2a, one gets the equation 

(e - 2a) (e + 2a) = 
with two possible solutions e = ib2a that give the differential equations 

{x^ - 4) y"{x) + 2{x + l) y'{x) - n{n - 3) y{x) = , (57) 

and 

(a;2 - 4) y"{x) + 2 (a; - 1) y {x) - n{n - 3) y{x) = . (58) 
They correspond to the densities 



respectively, hence the orthogonal polynomials are multiples of the Jacobi polynomials 

Pi^/^'-^/^)(r./2)andp(-^/^'^/^)(:r/2). 

Finally, in the third of the above cases, i.e. for d = 3a, we get again e = and 

{x'^ -4)y"{x)+^xy'{x)-n{n-4)y{x) = (59) 
corresponding to the density 

p(x) = — ''(^) ^ = v4 — a;^ . 

a{x) 

The corresponding orthogonal polynomials are multiples of translated Chebyshev polynomials 
of the second kind 

Pnix)= kn Snix) = ^ Snix) = ^ Unix/2) (n > 0) (60) 

n + 1 n + 1 
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(see e.g. [1], Table 22.2, and (22.5.13); Sn{x) are monic, see also Table 22.8), hence 

Pn{x) = pn-l{x) = — Un-l{x/2) (n > 1) . 

n 

We see that the recurrence equation (45) has four different classical orthogonal polynomial 
solutions! 

Example 2 As a second example, we consider the recurrence equation 

Pn+2{x) -{x-n - l)pn+l{x) + a {u + lfpn{x) = (61) 

depending on the parameter a G ffi. Here obviously the question arises whether or not there are 
any instances of this parameter for which there are classical orthogonal polynomial solutions. 
In step 6 of Algorithm 1 we therefore solve also for this unknown parameter. This gives a 
slightly more complicated nonlinear system, with the unique solution 



b=2c,c = c,d = —4 c, e = 0, a = 0, a 



Hence the only possible value for a with classical orthogonal polynomial solutions is a = 1/4, 
in which case one gets the differential equation 



+ ^) Pni^) - 22:^(2;) -2npn{x) = 



with density 

in the interval [—1/2, 00], corresponding to shifted Laguerre polynomials. 

4 Application: The Legendre Addition Theorem 

As an application of Algorithm 1 in this section we show how the particular case 

Pn{x^ + {l-X^)cOSe)=Pn{xf + 2j2T—W^Pnixf COS kO . (62) 

k=i y^^ + ^i- 

of the Legendre addition theorem ([14], 5.4.4, p. 239, see also [22], [11]) 

Pnixy + Vl-x^^Jl-y^ COS e) = Pn{x) P^iv) + 2 Pilix) P^{y) cos ke 



k-- 



can be deduced by linear algebra techniques. Note that (62) played an essential role in 
Weinstein's proof of the Bieberbach conjecture [22]. Here Pn{x) = Pji'''^\x) are the Legendre 
polynomials, and P^ix) are called the associated Legendre functions. Our goal will be to 
identify these functions. In our deduction, we partially follow [6], see also the first author's 
review [9]. For the given purpose, we write 

Pnix"^ + (1 - x"^) COS B) = B°{x) +2^2 r^— TTT (1 - a;^)^ ^nix) COS kO (63) 

k=l + 
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with still unknown functions B^{x). Multiplying by z", and summing for n = 0, 1, . . . yields 
the generating function of the Legendre polynomials, hence ([1], (22.9.12)) 

^1-Z{2X^ + {1-X^){W+IH)+Z^ = S £n " ^^^^ 

where we put w = e*^ and D~^{x) = B^{x). In the sequel we consider this equation as a 
formal Laurent series expansion w.r.t. the variables w and z. The functions BI^{x) can be 
iteratively calculated by series approximations of the left hand function (e.g., using Maple), 
and it turns out that, for < A; < n < 10, for example, these form polynomials that are 
squares of another system of polynomials 

Dt{xf = B'jx). (65) 

We normalize D^{x) such that the highest coefficient has sign (— l)'^ (to be consistent with 
the definitions given in ([1], § 8)). 

Now, we would like to find a three term recurrence equation w.r.t. n valid for the polynomials 

For this purpose, we "guess" that 

{ak + bn + c) D^^^i^) + {dk + en + f) D^^^{x) + {gk + hn + i) D^{x) = 

with unknowns o.,b,c,d,e, f,g,h,i. Substituting the given values D^{x) {0 < k < n < 10) 
into this proposed recurrence equation yields a linear system which turns out to be consistent 
(although we have (2) = 36 equations, but only 9 unknowns), with the unique solution 

in-k + 2) D^^^ix) - (2n + 3) x D^^^{x) + (n + k + I) D^{x) = . (66) 

Currently this recurrence equation is not yet proved, but this will be done soon. Assume for 
the moment that E^{x) are solutions of (66). Then by another application of linear algebra 
(see e.g. [20], [12]), this recurrence equation can be "squared", i.e. it is possible to calculate 
the recurrence equation of third order valid for E^{x)'^. This step can be accomplished, e.g., 
by the procedure rec*rec of the gfun packacke ([20], see also [10]) with Maple, and results 
in the recurrence equation 

{2n + 5){k + n + 2){k + n + l f S^{x) 

- {2n + 3){k + n + 2){k'^ -n^ + 4x'^ - 4n + 16 a;^ n + 15 a;^ - 4 ) 5^+i(a;) 

- {2n + 5){-2 + k-n){k'^ -n^ + 4:x'^n'^ -4:n + lQx'^n + 15x'^ -4)5^+2(^) 

+ {2n + 3) {-2 + k - n) {-3 + k - nf S^+six) = (67) 

for the squares S^{x) = E^{x)'^. 

If we are able to prove that this recurrence equation is valid for our unknown functions B^{x), 
then, by an a posteriori argument, we have deduced (65), since we have luckily found the 
"square root" recurrence equation (66) of (67). 

Next we show how it can be discovered independently that B!^{x) satisfy (67). We can rewrite 
(64) as 
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with 

n'^lz ■= ^ i 

' ^ ■ ^l-z{2x^ + {l-x^){w + llw))+z^ z^w^ 

where CTz,wG^{z^w) denotes the constant term of the double Laurent series G^{z^w). To 
obtain a recurrence equation for F^{x), we try to find polynomials Pi) ^25^3 in the variables 
n, A;, and a;, and polynomials Gi, and G2, both of degree 2 in both z and lo, such that 

Substituting Gi and G2 generically, and dividing by G^, a polynomial identity is derived, and 
by equating coefficients w.r.t. z and u;, we get a linear system, again. Solving this system 
results in the identity 

-{n + ?,){k + n + 2){2 + n-k) G^{z, w) 

+ ( n + 2 ) ( 4 rr^ - + 22 n - 6 n - 9 + 30 a;^ + A;^ ) G^+i {z, w) 

- (4n2 a;^ -n^ + 18 na;^ - 4n + 20a;2 + A;^ - 4 ) ( n + 3 ) G^+2(^: 
+ (n + 2)(n + 3 + A;)(n + 3-A;)G^+3(;2,i(;) = 

. d f G,{z,w)G^Jz,w) \ ^ d ni + z)G2{z,w)G^Jz,w) \ ^^^^ 



dz \ z^ w I dw 



z" w 



for certain polynomials Gi{z,w) and G2{z,w) which are reproduced in the appendix. 
Since any formal Laurent series / (z) satisfies 

applying CT^^^u to the identity (68) yields the three term recurrence equation 

-( n + 3 ) ( A; + n + 2 ) ( 2 + n - A; ) F^(a;) 

+ ( n + 2 ) ( 4 a;^ - + 22 n a;^ - 6 n - 9 + 30 a;^ + A;^ ) F^+i (a;) 
- (4n2 a;^ -71^ + 18na;2 - 4n + 20a;2 + A;^ - 4) (n + 3) F^+2(a;) 
+ (n + 2)(n + 3 + A;)(n + 3-A;)F^+3(a;) =0 

for F^{x). Using the recurrence equation 

(1 + A; + n) a„ + (-1 + k - n) a„+i = 

which is valid for the factor (n + A;)!/(n — A;)!, one can finally show (e.g., with rec*rec again) 
that B^{x) satisfy the recurrence equation (67), as announced. 

Now we know that B!^{x) are the squares of the polynomial system D^{x) defined by (66). 
But who are these polynomials? A direct application of Algorithm 1 is not possible since our 
polynomials D^{x) live for n > A; rather than for n > 0. Hence, we deal with 

T^{x) ■.= Dl^,{x) 

for which, by (66), we get the recurrence equation 

(n + 2) T^+2{x) - {2n + 2k + 2,)xTI^^^{x) + (n + 2A; + 1) tI^{x) = . 
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Now an application of Algorithm 1 shows that the differential equation 

{x'^-l) Tf{x) +2{l + k)xT!^'{x)-n{n-2k-Z) T^{x) = 
is valid, and that = ^"^^^"'"^ ■ Hence one has the density 

p{x) = {l-x^)\ 
i.e., we have multiples of Jacobi-Gegenbauer polynomials 

n\ 

with some function K{k) not depending on n, and 

1 " n 

Q(",/3) (^) = 1 f " + + f ) + 1)^ _ 1)"-^ 

^2n + a + /?-j^V fc ) \ n- k ) ^ ^ ' ^ > 

denoting the monic Jacobi polynomials (see e.g. [1], (22.3.1)). 

To determine k(A;), we compare the coefficients of cos" in (63). Since ([1], (22.3)) 

(2fl)! ^^jj ^ 1\ „„j rp I \ on— 1 ™n I rr.n—^\ 



Pn{x) = ^^x^-rO{x''-') and T„(a;) =2"-^a;" + 0(a; 



we get 



P„(a;2 + (l-a;2)cos0) = frc^ + (1 - a;^) cos + 0(cos"-i(^)) 

(2n)! 



2"n!2 



(1 -a;^)" cos"0 + O(cos"-^(0)) 



^ {I - x^Y D'f^ixf cosne 



{2n)\ 

^ (1 - x'^Y Ki^f 2*^"^ cos" 6 + O(cos"-i(0)) 



(2n)! 

because cos nO = Tn{cos6). Hence, we arrive at 

j^n( .2 ^ (2")!^ (l-x'^r 

and, setting n = k, we see that we must choose K{k) = (— l)'^ according to the nor- 
malization with highest coefficient sign (—1)*^ of D^{x). We have finally determined the 
representation 



and therefore 



^ (-1)' (1 - .-r-' Qt%) . 



This finishes the proposed identification of P^ix) 
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5 Classical Orthogonal Polynomials of a Discrete Variable 

In this section, we give similar results for classical orthogonal polynomials of a discrete variable, 
see Chapter 2 of [19]. These are given by a difference equation 

a{x) AVy{x) + t{x) Ay{x) + A„ y{x) = , (69) 

where 

Ay(a;) = y{x + 1) - yix) and Vy(a;) = y{x) - y{x - 1) 

denote the forward and backward difference operators, respectively. The classical discrete 
orthogonal polynomials are the polynomial solutions of the difference equation (69). 
These polynomials can be classified similarly as in the continuous case according to the func- 
tions (j{x) and t{x); up to linear transformations the classical discrete orthogonal polynomials 
are classified acccording to Table 2 (see [19], Chapter 2, and § 6). We added the other pos- 
sible solutions to those given in ([19], Chapter 2). In particular, case (2a) corresponds to the 
non-orthogonal solution x"' in Table 1. Similarly as for the powers 

-^x'' = nx''-^ , 
ax 

the falling factorials x— := x{x — 1) ■ ■ ■ {x — n + 1) satisfy 

Ax^ = nx^. 



It turns out that they are connected with the Charlier polynomials by the limiting process 

J = {-ir{x-n + l)n=x^ 



lim(-l)"/i"cl'')(a;) = lim(a;-n + l)„ ifJ 



where we used the hypergeometric representation given in ([19], (2.7.9)). 
Note, however, that other than in the differential equation case the above type of difference 
equation is not invariant under general linear transformations, but only under integer shifts. 
We will have to take this under consideration. 

The classical discrete orthogonal polynomial systems correspond to a discrete weight function 
p{x), with ([19], Equation (2.4.1)) 

+ 1) ^ (j{x) + t{x) ^^^^ 
p{x) a{x + 1) 

given by the Pearson type difference equation ([19], Equation (2.1.17)) 

A(a{x) pix)) = t{x) p{x) . (71) 

As in the continuous case, by Lemma 1 (the measure p is discrete now) they satisfy a recurrence 
equation (4) 

'Pn+l{x) = {AnX + Bn)Pn{x) - CnPn-l{x) ■ 

Multiplying (69) by p{x), by (71) the difference equation takes the selfadjoint form 
([19], Equation (2.1.18)) 



A (cj{x) p{x) Vy{x)) +Xny{x)=0 
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1. a{x) = 1, t{x) =ax + f3 =^ Pn{x) = Kn ' (x), see § 6, 

2a. (j[x) = X, (j[x) + t{x) = =^ Pn{x) = X— := x{x — 1) ■ ■ ■ {x — n + 1), 

2b. cr{x) = a;, cr{x) + t{x) = /j, {/j, ^ 0) =^ Pn{x) = c^\x), the Charlier polynomials, 

3 a{x) = x, a{x) + T{x) = jj, {'J + x) =^ Pn{x) = mn'^\x), the Meixner polynomials, 

4 a{x) = X, a{x)+T{x) = {N—x) =^ Pn{x) = kn\x,N), the Krawchouk polynomials, 

5 a{x) = x{N + a-x), a{x) +t{x) = {x + P+1){N -1 - x) ^ pn{x) = hl^'^\x, N), the 
Hahn polynomials (of first kind), 

6 a{x) = x{x + ij), (j{x) + t{x) = {iy + N-l-x){N-l-x) ^ pn{x) = h\^'''\x,N), the 
Hahn polynomials (of second kind). 

Table 2: Normal Forms of Discrete Polynomials 



I' 



from which one can deduce (as in Lemma 2) that Pn{x) satisfy a difference rule (see also [19], 
(2.2.10)) 

a{x) ^Pn{x) = anPn+l{x) + PnPn{x) + -fnPn-l{x) ■ (72) 

Analogous to the continuous case, one has a Rodrigues formula ([19], (2.2.8)) 

Pn{x) = (^pix) n ^(^ - (73) 

for some En independent of x, and the orthogonality relations read as ([19], (2.3.1)) 

(0 if m 7^ n 

/i„ / if m = n 
for some discrete set X C M. 

To deduce formulas similar to those in Theorem 1, we substitute (25) 

Pni^x^ = X -\- X -\- kj^ X -\- . . . 

in the difference equations for Pn{x) and for p„_|_i(a;), in the recurrence equation and in the 
difference rule. Equating the three highest coefficients in any of these four equations yields 
twelve nonlinear equations in the twelve unknowns 

■^m Bji, Cjii OLm Pm Im -^n+l: k^, k^, k^^i ■ 

In particular the highest coefficient of the difference equation yields 

Xn = —{an{n — 1) + dn) 

again. We can assume that A„ / for n G N, hence a{n — 1) + d / for n G N since otherwise 
no orthogonal polynomial solutions can exist; in particular, d ^ 0. 

By a Grobner basis calculation (invoked by the solve command of Maple or REDUCE) it 
turns out that there is a unique solution of the above nonlinear system, given by 
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Theorem 2 Let Pn{x) = knx"' + ... (n G No) be a family of polynomial solutions of the 
system of difference equations (69) that are orthogonal with respect to a discrete weight 
function p{x). Then the difference rule (72) 

a{x) Vpn{x) = anPn+l{x) + (3nPn{x) + 'JnPn-lix) 

is valid with 

an = an- 



(in 



In 



n {a{n — 1) + d) (2(n — 1) a {a{n — 1) + a + d) + ad + 2ae — hd) 
{2an + d) (2a(n - I) + d) 

n {a{n - I) + d) {a{n - 2) + d) 



-1 (a(2n - 3) + rf) (a(2n - I) + d) (2 a{n - 1) + df 



■ {^{n-l)[a{n-l)+d)[a^[n-lf + a[n-l)d+Aac+2ae-bd-b^)+ae'^-bde + cd^^ , 
and the recurrence equation (4) 

Pn+l{x) = {AnX + Bn)Pn{x) - CnPn-l{x) 

is valid with 

An = ^^, (74) 



_kn+i {n-l){d + 2b){a{n-l)+a + d)-2ae + d'^ + de + 2hd 
''~~kn {2an + d) {2a{n-l) + d) ^ ' 

and 

Cn = -^ , □ (76) 

kn a{n — 1) + d 

Note that with an immense effort one can deduce Theorem 2 also by hand calculations, using 
a similar technique as in our proof of Theorem 1 . 

Theorem 2 is also valid for case (1) of Table 2, see § 6; the recurrence equation part is also valid 
in case (2a) of Table 2, and generates the recurrence equation pn+2{x) — {x — n — 1) pn+i{x) = 
for the falling factorial Pn{x) = x—. 

As a byproduct of Theorem 2, we get for the ratio 

Corollary 6 Let Pn{x) = knx'^ + k'^ x'^^^ + . . . (n G No ) be a family of orthogonal polynomial 
solutions of the system of difference equations (69). Then the relation 

A;^+i_n + l {dn + 2hn + 2e){2a{n-l) + d) kn+i 
k'n ~ n {d + 2an){2b{n-l)+dn + 2e-d) kn 

is valid. □ 

As in the continuous case, the Rodrigues constant En, given by (73), can be determined. 
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Corollary 7 Let Pn{x) = A;„ + . . . (n G No) be a family of orthogonal polynomial solutions 
of the system of difference equations (69) corresponding to a discrete weight function p{x). 
Then pn{x) have a Rodrigues representation (73), En being a hypergeometric term, satifying 
the relation 

En+i {a{n -l)+d) kn+i 



En {a{2n - 1) + d) d kn 
Proof: In ([19], (2.2.10)) it was shown that 



(77) 



a{x) Vpn{x) = ^" (Tn{x)pn{x) - Pn+A (78) 

where T„ (a;) = t{x + n) + a{x + n) — a{x) ([19], (2.1.15)). An application of a method similar 
to that used in Theorem 2 yields the coefficients for this type of difference rule, and hence 
determines En+i/En according to (77). □ 

Also analogously we get 

Corollary 8 Let Pn{x) = kn x""' + . . . (n G No) be a family of orthogonal polynomial solutions 
of the system of difference equations (69) corresponding to a discrete weight function p{x). 
Then the relation 

= {d + an — a){—a^ n'^ — 4dcna + dbe + db^ n — d^ c — ae^ — 4cn^ 



2nade-2n^ d + n^ adb + bnd^ -n^ ad^ +r?b'^ a-2a^ er?){n+l) 

fk ^ ^ 

l({d + 2an + a){d + 2an- a){d + 2anfY "^^ 

is valid. 

Proof: In ([19], (2.5.6)) it was proved that 

Cn kn—A hfi. 



hence 

^71+1 Cn+l ^71+1 



hn ^71+1 f^n 

An application of Theorem 2 yields the result. □ 

Similarly as in the continuous case, Theorem 2 can be used to generate an algorithm to 
test whether or not a given holonomic recurrence equation has classical discrete orthogonal 
polynomial solutions. 

Algorithm 2 This algorithm decides whether a given holonomic three-term recurrence equa- 
tion has classical discrete orthogonal polynomial solutions, and returns their data if applicable. 

1. Input: a holonomic three-term recurrence equation 

qn{x)Pn+2{x) + rn{x) Pn+i{x) + Sn{x) Pn{x) = (g„(a;), r„(a;), s„(a;) G Q[n,x]) . 

2. Shift: Shift by max{n G No |n is zero of either g„_i(x) or Sn{x)} + 1 if necessary. 
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3. Rewriting: Rewrite the recurrence equation in the form 

Pn+i{x) = tn{x)pn{x) + Unix) pn-i{x) (i„ (x) , it„ (x) GQ(n,a;)) . 

If either tn{x) is not a polynomial of degree one in x or Un{x) is not constant with respect 
to X, return "no classical discrete orthogonal polynomial solution exists"; 
exit. 

4. Linear Transformation: Rewrite the recurrence equation by the linear transformation 
X with (yet) unknowns / and g. 

5. Standardization: Given now A^^B^ and C„ by 

Pn+l{x) )Pn{x) - CnPn-l{x) (^n, -Bn, C^n G Q(n), / 0) , 

define 



n 



according to (74). 
6. Make monic: Set 

Br,:=^e Q{n) and (?„ := G Q(n) 

and bring them in lowest terms. If the degree of either the numerator or the 
denominator of -B^j is larger than 2, if the degree of the numerator of C„ is larger than 6, 
or if the degree of the denominator of Cn is larger than 4, then return 
"no classical discrete orthogonal polynomial solution exists"; exit. 



7. Polynomial Identities: Set 

according to (75), and 



Bn — , Bn 



— ~, 



according to (76), in terms of the unknowns a, 6, c, c?, e, / and g. Multiply these identities 
by their common denominators, and bring them therefore in polynomial form. 

8. Equating Coefficients: Equate the coefficients of the powers of n in the two resulting 

equations. This results in a nonlinear system in the unknowns a, 6, c, d, e, / and g. Solve 
this system by Grobner bases methods. If the system has no solution, then return 
"no classical discrete orthogonal polynomial solution exists"; exit. 

9. Output: Return the classical orthogonal polynomial solutions of the difference equa- 
tions (69) given by the solution vectors (a, 6, c, e, /, g) of the last step, according to the 
classification given in Table 2, together with the information about the standardization 
given by (74). This information includes the necessary linear transformation jx + g, as 
well as the discrete weight function p[x) given by (70) 

p{x + 1) _ c!{x) + t{x) 
p{x) o{x + 1) 
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Proof: The proof is an obvious modification of Algorithm 1. The only difference is that 
we have to take a possible linear transformation fx + g into consideration since the differ- 
ence equation (69) is not invariant under those transformations. This leads to step 3 of the 
algorithm. □ 

Note that an application of Algorithm 2 to the recurrence equation pn+2{x) — [x — n — 1) 
Pnj^i{x) = which is valid for the falling factorial p„(j;) = a;—, generates the difference equation 
xAVpn{x) — xApn{x) + npn{x) = of Table 2:2a. 

Example 3 We consider again the recurrence equation (61) 

Pn+2{x) - {x-n - l)pn+i{x) + a {n + lfpn{x) = 

depending on the parameter a G IK. This time, we are interested in classical discrete orthogonal 
polynomial solutions. 

According to step 3 of Algorithm 2, we rewrite (61) using the linear transformation x i-^ 
with yet unknowns / and g. Step 4 yields the standardization 

kn+l _ 

f 

In step 7, we solve the resulting nonlinear system for the variables {a, 6, c, d, e, /, g, a}, resulting 
in 

/ n A A b{-e + d+b) d+2b e b{d+b) \ 

This is a rational representation of the solution. Since we assume a to be arbitrary, we solve 
the last equation for b. This yields 

2 V ^/^^y 

which cannot be represented without radicals. Substituting this into (79) yields the solution 

f ^ , (i / 1 \ 4:ae — e — 2ad el. 1 el 
<a = 0, 0= 1± , , c = ; r — lb , f = =F / =1 Q = f , 

I ' 2V VT^j 2(i-4a) ^yr^'^ dr 

d and e being arbitrary. It turns out that for 1 — 4 a > this corresponds to Meixner or 
Krawchouk polynomials. 

Example 4 Here we want to discuss the possibility that a given recurrence equation might 
have several classical discrete orthogonal solutions. Whereas the recurrence equation of the 
Hahn polynomials hl^'^\x, N) has (besides several linear transformations) only this single 
classical discrete orthogonal solution, the case (3 = —a results in two essentially different 
solutions. 

Here one has the recurrence equation 

{n + 2 + a) {2 + n) {2n + 2) {n - N + l)pn+2{x) + (3 + 2n) 

{-Qua -211"^ a - An"^ X - I2nx + 2n^ N + 6nN + 4N - 4a - 8x) (x) 
- (1 + n) (n + 1 - a) (2n + 4) (n + TV + 2) p„(a;) = . 
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An application of Algorithm 2 shows that this recurrence equation corresponds to the two 
different difference equations 

X {-X + l-a + N)^S.Vpn{x) + {-2x + N + aN) Apn{x) ) + n (n - 3)pn{x) = 

and 

{x + a) {-X + 1 + N) A\7p„{x) -{2x-N + 2a + aN) Apn{x) +n{n- S)p„{x) = . 



6 A New Polynomial System 

In this section, we would like to present a new polynomial system satisfying the difference 
equation (69), but not coming from a discrete weight function, see Table 2:1. We set a{x) = 1, 
and t{x) = ax + (3 {a ^ 0). This yields the polynomials 



1 A f —n r + 



a 



-ir-2Fo 



a 



where we have used kn = a". Note that the two different hypergeometric representations 
come from ([19], (2.7.7)-(2.7.8)): they convert into one another by changing the direction of 
summation. The functions K^'^\x) are not covered in ([19], Chapter 2) since 



p{x) 



l + ax + j5, i.e. p{x) = { ^^ ^ \ 

\ a J. 



does not correspond to a valid weight function p{x) over a suitable discrete real set. 

The deductions in [19] show that (77) is valid, and we find that the Rodrigues constant is 

En = 1, hence (73) 

p{x) ((l + /3)/a), V \ a Jj 

Applying Zeilberger's algorithm ([24], see also [7]) to one of the hypergeometric representations 
yields the recurrence equation 

{x)-{na + xa + a + P) if^^f {x)-{l+n)a K^^^f^^ {x) = (80) 

with the initial functions 

K^^f\x) = {), K^o"'^\x) = l. 
The next three polynomials are 

k['''^\x) =ax + P , K^^''^\x) = x^Q^ + a {2 ^ + a) x + 0^ + a + a ^ , 

[x) = a^x^+2, a^{a+l3)x'^+a(& a (3+2 +3 (3'^ + 3a)x+(3^+2 0^/3+3 a /3^+3 a (3+2 c? . 
According to (78), the difference rules 

VK("'^)(2;) = -(ax + an + /3)K("'^)(x)+ifi^f (a;) =ani^i'!.f (a;) (81) 
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are valid which are consistent with Theorem 2. Another appHcation of Zeilberger's algorithm 
yields furthermore the recurrence equation 

{I + (3 + xa + a) K^'^'^^x + 2) - {(3 + x a + na + 2 + a)K^'''^^x + I) + K^'^'^^x) = 

with respect to x which also can be obtained combining (80) and (81). 

Note that the given polynomial system completes the classification of the polynomial solutions 
of (69) given in ([19], Chapter 2), resulting in Table 2. Since the given system satisfies a 
difference rule which is compatible with Theorem 2, Algorithm 2 can also recognize these 
solutions. 



7 Coefficient Recurrence Equation 

By equating the coefficients of the powers of x in their differential equation (17) one gets for 
the continuous classical orthogonal polynomials 

n 

Pn{x) = J2a^^^x' (82) 

the three term recurrence equation (A; = 0, . . . , n — 1) 

{n-k){a{n + k-l) + d) af^ = {k + 1) [{k h + e) ag, + (A; + 2) cag^) , ^ 0) 

which was given in ([13], Eq. (3)). 

In this section, we discuss the same question in the discrete case. Since the operators A and V 
behave nicely with falling factorials only, the solution is a little more difficult in this setting. 
Throughout this section we assume an expansion (82), i.e., we use the renamings kn = a'n\ 

Let now Pn{x) be a classical discrete system satisfying the difference equation (69) 

a{x) AVpn{x) + t{x) Apn{x) + XnPn{x) = 

with a{x) = ax^ + hx + c^ t{x) = dx + e and A„ = — (an(n — 1) + dn). Note that the difference 
equation is equivalent to 

{a{x) + t{x)^ Apn{x) = a{x) Vp„(x) - A„p„(a;) , (83) 

hence the use of the double difference operator A V can be omitted. 
Taylor's theorem yields the expansions 

k=l ' k=l 

and by Leibniz's formula we receive the representations 

n—ln—k , n—ln—k 

Ap.(:.) = E E ( t ' ) "S^- vp„(.) = E E ( t ' ) (-1)^'^^ 

k^O j=l k^O j=l 

for Apn{x) and Vpn{x). 
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Substituting these representations in (83) and equating the coefficients of the powers of x 
results finally in the recurrence equations 

(2a(n -l)+d) a^:\ = n [e + [h +'^){n - 1)) , (84) 

and for < A; < n — 2 

{n-k){a{n + k-l)+ d) af'^ = ((26 + d)k + 2e) aJ^'J^ + (85) 

Note that from (84) and (85) one can deduce the discrete equivalents of (29) and (33). 

As an example, we apply these formulas to the falling factorials Pn{x) = x- which are solutions 

of the difference equation (Table 2:2a) 

x/S.Vpn{x) - xApn{x) + npn{x) = . 

Here we have a = c = e = 0, b = —d = 1. Hence for the coefficients s{n, k) of 

n 

x-=^ s{n, k) x^ 
that are called the Stirling numbers of the Brst kind, we get the relation 

n—k 

{n-k)s{n,k) = Y.{-iy{':ll)s{n,k+j) 
which is due to Lagrange (see e.g. [5], p. 215, Theorem C). 

Appendix 

Using higher coefficients k'^,k'^', ... of the polynomials 

y{x) = pn{x) =knx'' + k'^ + k'J^ a;"-^ + k'^^ x""-^ + C ^""^ + • • • (n G Nq, A;„ / 0) 
that are solutions of (1) 

a{x) y"{x) + t{x) y'{x) + A„ y{x) = , 

or (69) 

(j{x) A\7y{x) + t{x) Ay{x) + A„ y{x) = , 

respectively, it is possible to use these as auxiliary variables, hence determining formulas for 
them. Note that to obtain the results of this article, it was crucial to make use of k" which 
was never done before. 

In this appendix, we collect results for the coefficients k'^,k'^ and k'^' in both the continuous 
and discrete cases. 



31 



In the continuous case, we have the following properties: 

^■"+1 _ ( iP + 21) e. n + 2 c n a - iP n - h e. + c d + ) 



and 



k'^ [rfilP -3lPn + 2ben + 2cna + cd + 2lP -2ca-3he + e^) 

{n + \ ) {2an — 2a + d) {d — "ia + 2an) k^+i 
{d — a + 2an) {d + 2an) {n — 1) kn ' 



= (n + l)(n^ -Sn'^b^ + en'^bca + 3n'^b'^e-enbca + 3nbcd + eneca 

-Gnb'^e + Snbe'^ + 2nb^ + 3ecd-2bcd-3be'^ -2cea + 2b'^e + e^) 
{d - 3a + 2an) {2an - Aa + d) j {^{ d + 2an) {d — a + 2an) {n — 2){n b 

+ 6 6ca + 3 6^ e-6 6^ + 6neca + 11 n 6^- 18 n6ca - 12 n 6^ e 
+ 3nbe^ + 3nbcd-Qb^ -hbcd-Qbe^ + llb'^ e + 3ecd-^cea 



+ 126ca + e^)] 



-Mi = (n + 1)(6 5^ cd+ 1171^6^ + 12 c^n da + nH^ -6 nH^ -36 n^ft^ca 

kr, 



"n 



+ I2nce^ a + 22nb^ e - l%nb^ + Anbe^ - b^ e + & b^ + 4.11^ b^ e 

- Mnbcea + 12nebcd - 8ce^ a - lAebcd + 24:nb'^ ca + 6 6^ cd - 6be^ 
+ 6 6^ cd + 11 6^ - 6 6^ e + 12 6ce a + 3 ^2 _ 6 da - 12 n + e'' 
-6nb'^ + 12c'^n'^a'^ +24:n'^bcea-Unb^ cd+12n^b^ ca){d-6a + 2an) 
{2an-4:a + d) ^(^{d + 2an){d- a + 2an){n-3){6n'^b'^ cd + 35n'^b'^ 

+ 120"^ nda + n^b'^ - lOn^ b'^ - 12r? 6^ ca + 12 n ce^ a + 70 n 6'^ e - 30n 6^ ^ 
+ Anbe^ -30n'^b^e + 671^ e'^b'^ + 4:n^b^e-92nbcea + 12nebcd 

- 20 ce^ a - 26 e 6cd + 132 n 6^00 + 6 e^cd- 10 + 26 6^0^ + 35 6^ 

- 50 6^ e + 80 6ce a + 3 - 18 da - 36 n + e'' - 50 n 6"^ + 12 

+ 24n26cea-26n62cd+ 12nH2ca + 24c2a2 - 7262ca + 246'')) . 

Even more complicated results are valid in the discrete case: 
k" 

= (n + l)(2a2 n^ + 12d6n2 - 6 + 3 d^ + 5 a dn^ + 12 6^ + 24e6n 

+ 12 aen - d^ n-7a dn + 24 can + 12 den + 4a^n - 12 d6n - 12 6^ n 

- 12 6e + 2ad-2 d^ + 12 cd + 12 e^)(2an-2a + d)(d-3a + 2an)/( 

(d + 2an)(d — a + 2an)(n — 1)(12 den + 5adn^ + 12 aen — 36 d6n 

+ 12 d6n^ + 24 e6n + 2 d^ + 12 + 24 6d + 22 a^n- 12 a^n^ - 12a e 

+ 3 d^n^ - 7d^n - 12 de + 2 a^n^ + 12 cd + 24can - 36 6^ n - 36 6e + 14a d 

kn+l 



+ 246^ - 12a2 - 24ca + I2n'^ 6^ - 17adn)) 
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and 



-III 



I^±i = (n + l)(20ea2„ + 46„4^2 _ -^g^^^a - ^ dn^ - 2en<f + Aen^ a'^ 



+ 2dn^ a? - %h(? AhnS - %er? + 12eda - 13 r? a-Adnc? 
+ 5 a + 10 d + 6 a n + 20 6 - 36 6^ d - 8 e + 12 6 d a n 

- 8ae^ + 16n6^ + 8nH^ +n^d^ - 2471^6^ - 26eand + 22en2ad 
+ 10 6n^ad- Mfon^d^ - 26 6n2da + 24rid62 - 24n6ae - 36n6de 

+ 24n2ea6 + 86(ie + Gn^d^e + 12nde2 + 12nd2c + 6n^d^6 + 24nae2 
+ 24neH + 24n^e6^ + 8 -8 cd^ -24 fee^ + 16 6^6 + 24 n^de 6 
+ 24n^ dca + I2n^ db^ + 24e cd - ASneb^ - 16bcd -24ndca + 48neca 
+ 24nbcd-48nbca + 48n'^bca-l6cea-d^n'^ -4ed'^ -2d^n) 
{d-3a + 2an){2an-4a + d) j {^{ d — a + 2an) {d + 2an) {n — 2){ 

48 ea^ n + 46n^ - 32 6n^a^ - 16 dn^a^ - 14e nd^ +4en^a^ + 2dn^ 

- 112 6 a^n + 50 6nd^ - 20 en^a^ + 60 e da - 28 d^n^ a - 56 dna^ + 5 d^n^ a 
+ 46 d + 47 d^ a n + 92 6 - 24 d^ a + 24 d - 72 nH^ d - 48 6 a d 

- 40 e + 94 6dan - 48 6^ - 32 a + 88 n 6^ + 8 6^ + d^ - 48 6^ 

+ 48 6a^ - 24 6d^ +48 6e a + 96 c6a - 12 e2d + 48 cda - 72 d6^ - 70 e and 
+ 22 en^ ad + 10 6n'' ad - 32 d^ - 56 da + 132nd6^ - 72n6 a e 
-84n6de + 24n2ea6 + 686de + 6n2d2e + 12nde2 + 12nd2c + 6n^d^6 
+ 24nae2 + 24716^6 + 2471^6 6^ + Be''' -20cd2 -48 662 + 88 6^6 
+ 24n2d66 + 24n2dea + 12 d62 + 246 cd- 96 n 6 6^ -40 6cd 

n6c( 



72 n d ca + 48 neca + 24nbcd — 144 n6ca + 48n2 6ca — 64cea 



4 d^ + 4 e d^ + 3 d^ n)) 



k"" 

We omit the very lengthy result for -g^ in the discrete case. Similar expressions can be 
obtained for higher coefficients. 

These results should indicate that the use of computer algebra is quite natural in the given 
context, and hand calculations seem not to be adequate for these computations. 
Finally, for the sake of completeness, we present the full polynomial certificates Gi {z, w) and 
G2{z, w) from § 4: 

G2{z,w) = ^(-18 u; + 15z + S0zx'^w-15zx'^w'^ -10 z'^w + 15 zw'^ -15 zx'^ + 4zw'^n^ 

O ^00 o 

— 15zw A; + 3l2;u; n — 16znkx +16znk — 4zn x +4zn k — 15zkx 

^ o o 

-Slznx —20zn x + 15 z k + 31 zn + 4zn + 20 zn —2zw — 16zw nk 

O O 000 0*^0 00 00 

+ 16zw nkx +4zw n kx — 4zw n x — 4zw n k + 20 zwn 

00 00 000 O 00 o 

+ 15 zw kx — 31 zw nx — 20 zw n x +62zwnx —18 2; wn —24z wn 

f)o 00 QO 00 

— 22 wn — 4wn — 36«;n + 40 2;wn x +8zwn' x —2zwn + 4k z wn 

— &k^ w - 4z^ wn^ + 10k^ z^ w - 4k^ zw - 4k'^ wn - 4zn^ kx^) jk 
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and 

Gi{z , w) = — -{5 z'^ + 3 z + 3 z x'^ w'^ — 3 z w'^ — 3 z x'^ — 32 k z w n + 6 k z'^ w — 4 z'^ k'^ n x'^ 

- 8k zwn^ - 30k zw + 2Ak zwn^ x"^ + 78k zvjx"^ - 6zv? k"^ + 9 zv? k 
-hzu? n-lOz"^ k^ x^ -8kwr? - 28kwn + lQz^ k^ + lbz^ k + 7 z^ n-^oz"^ x^ 
+ 2z'^n'^ + l^z'^nk-l^z^nkx'^ -Az'^r^kx'^ +4:z'^n^k+4:z'^k^n-lbz'^kx'^ 
-1 z^ nx^ -2 z^ x^ - A zk^ nx^ - 12 znkx^ + 12 znk + A zr? k + Azk^ n 

— 9 z kx^ — hznx^ — 2zr?x^ — Qz k^ + Q zk^ + 9zk + 5zn + 2zn^ 

+ 10 z'^ w"^ k^x^ + 4: z^ w"^ k'^ nx'^ + I6z'^ w'^ nk - I6z'^ w'^ nkx'^ - A 'u? k x^ 

- 10 xi?k^ + h z^ V? x^ -2 z^ ui^ + 15 z^ ui^ k-7 z^ui^n + A z^ r? k 

- ^z^ k^ n-\hz^ kx^ + 1 z^ nx^ + 2 z^ w'^n^x'^ + 6z w"^ k^ x^ 

+ ^zw^ k^ nx^ — h z^ + 12 zw^ nk — 12 zw^ nk x^ — ^zw^ kx^ + ^zw^ k 

— Azw"^ n — 2 z 'iv'^ — 9 z ii? k + h zw"^ nx^ + 2z w"^ x^ — 32kz'^wn x^ 

00 O O 00\/ 

— 8k z wn X — 30 k z w x — 24:kw + 4kz wn + 88 kzwnx — 4zn kx ) k 
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